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ABSTRACT 

We report the discovery and initial characterization of Qatar-2b, a hot Jupiter transiting a V = 13.3 mag K 
dwarf in a circular orbit with a short period, P\, = 1 .34 days. The mass and radius of Qatar-2b are Mp = 2.49 Mj 
and Rp = 1.14R], respectively. Radial-velocity monitoring of Qatar-2 over a span of 153 days revealed the 
presence of a second companion in an outer orbit. The Systemic Console yielded plausible orbits for the outer 
companion, with periods on the order of a year and a companion mass of at least several Mj . Thus Qatar-2 joins 
the short but growing list of systems with a transiting hot Jupiter and an outer companion with a much longer 
period. This system architecture is in sharp contrast to that found by Kepler for multi-transiting systems, which 
are dominated by objects smaller than Neptune, usually with tightly spaced orbits that must be nearly coplanar. 

Subject headings: planetary systems — stars: individual (Qatar-2, 2MASS 13503740-0648145) — techniques: 
spectroscopic 



1. INTRODUCTION 

Wide-angle ground-based photometric surveys, such as 
WASP and HATNet, have been effective for the identification 
of close-in exoplanets that transit their host stars. Although 
most of the systems found by these surveys are fainter than the 
stars targeted by Doppler surveys, they are still bright enough 
to allow confirmation and characterization of their transiting 
planets using follow-up photometric and spectroscopic obser- 
vations. 

Most of the more than 150 confirmed transiting planets 
are hot Jupiters. Of those found by ground-based surveys, 
only 12 are smal ler than Saturn, i.e. less than 9.4 Earth radii 
([Schneider 2011, as of September 2011). In contrast, most 
of the more than 1000 candidates identified by Kepler are 
smaller than Neptu ne, i.e. less than 3.8 Earth radii (e.g., see 
lLatham et al.11201 Tb . However, it will not be possible to de- 
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termine spectroscopic orbits for the vast majority of Kepler's 
small candidates, because the required velocity precision is 
beyond the reach of present capabil ities with instru ments such 
as HIRES on Keck 1 (e.g., see lBatalha e"tal]l201 11) . 

Part of the problem is that most of the Kepler candidates 
are 14th magnitude or fainter, a result of the fact that Kepler 
monitors a region that is only 0.25% of the sky. To find all 
of the nearest and brightest transiting planets, we now need to 
extend our photometric surveys to cover the entire sky. The 
prospects are good that such surveys will eventually be pur- 
sued from space, with missions such as TESS (an Explorer se- 
lected for Phase A by NASA) and PLATO (proposed to ES A), 
but the earliest any of those spacecraft could be launched is 
2016. 

In the meantime, at least three teams are working on 
ground-based photometric surveys designed to find the best 
systems for follow-up studies. An obvious strategy is to 
cover as much of the sky as possible (idealy all), so that no 
prime systems are missed. A complementary strategy is to 
target smaller stars, where smaller planets are easier to de- 
tect and characterize. The MEarth project has adopted both 
strategies by focusing on the coolest and smallest M dwarfs, 
targeting a few thousand of the nearest and brightest exam- 
ples all over the sky. This appro ach has already yielded GJ 
1214b (Charbo nneau et all 120091) . which lies in the "Super- 
Earth" transition region between Neptune and the Earth and 
has attracted a lot of interest because its favorable contrast 
with its host star allows studies of its atmosphere (e.g., see 
iBean et all 1201 11) . The other two efforts, the extensions of 
HATNet to HAT South and of WASP to the Qatar Exoplanet 
Survey (QES), have chosen to improve the capability of the 
more traditional strategy of a magnitude-limited wide-angle 
survey. Both HAT South and QES have implemented cameras 
with larger apertures, thus providing better photometric pre- 
cision and lower rates of contamination by faint stars close to 
the targets (enabled by the more favorable pixel scale result- 
ing from the longer focal lengths of the larger cameras). An 
important trade-off of this approach is that it requires more 
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detector pixels to cover the same area on the sky. 

A practical implication of finding smaller planets is that the 
follow-up Doppler spectroscopy needed to determine orbits 
and planetary masses will demand access to the most capable 
facilities for very precise radial velocities, such as HARPS on 
the European Southern Observatory's 3.6-m telescope on la 
Silla in the South, and HARPS-N, scheduled to come into op- 
eration in 2012 on the Telescopio Nazionale Galileo (TNG) 
3.6-m telescope operated for the Italian Institute of Astro- 
physics (INAF) on La Palma in the North. 

In this paper we report the discovery and initial charac- 
terization of Qatar-2b, the second transiting planet from the 
Qatar Exoplanet S urvey to be confirmed (for the first, see 
lAlsubai et all 1201 lb . Qatar-2b is especially interesting, be- 
cause our radial-velocity monitoring shows that there is a sec- 
ond companion in the system with an orbital period of about 
a year and a mass of at least several Mj. 

2. QES DISCOVERY PHOTOMETRY 

The Qatar Exoplanet Survey (QES) is a wide-field photo- 
metric survey for transiting planets. The initial 5-cameraCCD 
imaging system is now deployed at a site in New Mexico and 
has been in operation for more than two years. It uses an array 
of five Canon lenses equiped with 4K x 4K CCDs to image 
an 1 1 x 1 1 ° field on the sky s imultaneously at two different 
pixel scales (Alsubai et al! l201 lb . The data are reduced at the 
University of St. Andrews using pipeli ne so f tware based on 
the image-subtraction algorithm of Bramich (2008), and the 
data products are archived at the University o f Leicester, us- 
ing th e same architecture as the WASP archive (Pollac co et al.l 
2006). 

The QES photometry for the V = 13.3 mag K dwarf 
that we now designate Qatar-2 (3U C 167-129863, 0:2000 = 
13 h 50 m 37?41, <5 2 ooo = -06°48'14"4 dZacharias et al.l 120101) ) 
revealed transit-like events, found by an automated search 
on the archive data using t he box least-squares algorithm of 
iKovacs. Zucker. & M azeh (20 021) as modified f or the Super- 
WASP project by ICollier Cameron et al.l (120061) . Systematic 
patterns of correlated noise were modeled and removed from 
the archive l ight curves using a combination of the SysRem 
algorithm of iTamuz. Mazeh. & Zucker (2005) an d the tren d 
filtering algorithm (TEA) of lKovacs. Bakos. & Noyesl (12005). 
3UC 167-129863 was screened and selected as a promising 
candidate, along with tw o other stars in the sam e QES field, 
using tests described by ICollier Cameron et al.l (120071) . de- 
signed to ensure that the depths and durations of observed 
transits are consistent with the expectations for objects of 
planetary dimensions transiting main-sequence stars. The dis- 
covery light curve, shown in FigureQ] indicated an ephemeris 
of E = 2455618.9092 ±0.0014 + JV x 1.337098 ± 0.000022 
BJD. The final photometric data for the QES light curves are 
provided in full in the electronic version of the journal. 

3 . FOLLOW-UP SPECTROSCOPY WITH TRES 

We monitored the spectrum of Qatar-2 using the Tilling- 
hast Reflector Echelle Spectrograph (TRES), mounted on the 
1.5m Tillinghast Reflector at the Fred L. Whipple Observa- 
tory on Mount Hopkins, Arizona. Over a period of 153 days 
from 201 1 January 18 to June 21 we obtained 44 spectra with 
a signal-to-noise ratio per resolution element (SNRe) of at 
least 20 in the continuum at the center of the order contain- 
ing the Mg b features (near 518 nm). We used the medium 
fiber, which yields a resolving power of R ~ 44,000, cor- 
responding to a resolution element with 6.8kms _1 FWHM. 



Th e spectra were extract ed using the procedures outlined 
by Buchha ve et al.l (I2010I) . The wavelength calibration was 
established using exposures of a Thorium-Argon hollow- 
cathode lamp illuminating the fiber, immediately before and 
after the stellar observations. A typical exposure time for a 
stellar observation was 45 minutes. 

3.1. Radial Velocities 

We derived radial velocities from the TRES spectra using 
the procedures for mu lti-order cross correlations described by 
iBuchhave et al.l d2010t) . We present here a brief summary of 
the technique. First, we rejected roughly half of the echelle 
orders: several orders with low SNR shortward of 446 nm; or- 
ders longward of 678 nm, which are contaminated by telluric 
absorption lines or show strong fringing in the CCD; and a 
few orders in between which are known to give poor velocity 
performance relative to the other orders (generally because of 
some other contamination such as interstellar Na D absorption 
or emission lines from the Earth's atmosphere). In all, we in- 
cluded 23 orders in our analysis. Each observed spectrum was 
cross correlated, order by order, against the corresponding or- 
ders from the strongest single observation, the one obtained 
on JD 2455646. The cross correlation functions (CCFs) from 
the individual orders were summed and fit with a Gaussian 
function to determine the radial velocity. We estimated the 
internal precision of the radial velocity for each observation 
by fitting a Gaussian to the CCFs of the individual orders and 
using the deviations from the global fit to calculate the stan- 
dard deviation of the mean. 

We adjusted the velocities to correct for small shifts in the 
zero point between observing runs, which mostly resulted 
from minor modifications to the hardware between some runs. 
We observed IAU radial-velocity standard stars every night 
and were able to establish the run-to-run shifts with an uncer- 
tainty on the order of 5 ms -1 . These adjustments have been 
applied to the individual velocities for Qatar-2 reported in Ta- 
ble [T] which are relative to the observation obtained on JD 
2455646. The absolute velocity of that observed template was 
derived using our observations of IAU radial-velocity stan- 
dard stars and correlations of the Mg b order against a library 
of synthetic templates, yielding a value of -23.8 ±0.1kms _1 
on a system where the velocity of our primary standard, HD 
182488, is defined to be -2 1.508 km s -1 . The error estimate 
for this absolute velocity is fairly large, because only a single 
echelle order was used in the comparison with HD 182488 by 
way of the synthetic library spectra, and the absolute velocity 
of HD 182488 itself has an uncertainty at the level of nearly 
0.1 kms -1 . 

3.2. Orbital Solution 

We ran a Markov Chain Monte Carlo (MCMC) analysis on 
the radial velocity data, starting with Gaussian constraints on 
the priors for period and epoch from the QES discovery pho- 
tometry. We found that a single Keplerian orbit matched the 
data poorly, with % 2 = 424.9 for 44 observations and 38 de- 
grees of freedom. Introducing two additional terms, for a 
linear and quadratic velocity drift in a Taylor expansion, re- 
sulted in a dramatic improvement to the fit, with \ 2 = 46.0 for 
36 degrees of freedom (reduced x 2 = 1.28). This model pro- 
vided convincing evidence that Qatar-2 is orbited by a tran- 
siting hot Jupiter, Qatar-2b, with period P = 1.34 days, plus 
a second companion, Qatar-2c, with a period of roughly a 
year. The velocity curves for the orbit of Qatar-2b (phased 
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FIG. 1 . — QES discovery light curve, phased on the transit ephemeris of Qatar-2b. The data shown here were obtained with Camera 403 in campaigns C2, C3 
andC4. 



to the period) and the acceleration due to Qatar-2c are plot- 
ted in the bottom two sections of Figure |2j together with the 
corresponding components of the observed velocities. The 
top panel shows the original velocities together with the com- 
bined velocity model. The panels immediately under the ve- 
locity curves show the residuals from t he fits. 

Next we used the Systemic Console dMeschiari et al.f 2009) 
to explore possible simultaneous fits for two Keplerian or- 
bits. A Levenberg-Marquardt solution gave a plausible fit and 
provided starting parameters for an MCMC analyis. When 
the eccentricities were allowed to be free parameters, sev- 
eral plausible solutions with similar parameters were found, 
all with periods a little shorter than a year and companion 
masses of several Mi. The key parameters for a typical solu- 
tion from the Systemic Console are reported in Table|2j along 
with the parameters from the quadratic solution. Note that 
the period of the inner orbit from the radial velocities alone, 
P RV = 1.337169 ±0.000076 is a close match to the period 
from the final global solution together with all the photom- 
etry, P = 1.3371 182 ±0.0000037 days. A reliable solution 
of the orbit of Qatar-2c will require additional observations, 
which are planned for the coming season. 

3.3. Bisector Analysis 

The radial velocity of a star is defined to be the velocity 
of the center of mass of the star along the line of sight. Ob- 
servationally, radial velocities are determined by measuring 
the Doppler shifts of spectral lines formed in the star's atmo- 
sphere. Distortions of line profiles due to phenomena such 
as dark spots on the surface of a rotating star or contamina- 
tion of the spectrum by another star with variable velocity can 
mimic radial-velocity variations in the target star. A common 



technique for detecting line profile varia tions involves mea- 
surements of lin e bisector spans (e.g., see lOueloz et alj|2001t 
iTorres et al. 2007). Any variation in the bisector spans that is 
correlated with variations in the radial velocities is a strong 
warning that the velocity variation is probably not due to the 
reflex motion induced by an orbiting planet. 

The variations in the bisector spans that we measured for 
the 44 TRES observations of Qatar-2 are reported in Table Q] 
and are plotted in the bottom panels of Figure [2] We were 
particularly interested to see if the bisector spans showed a 
correlation with the quadratic drift attributed to Qatar-2c, as 
might be expected if the outer companion contributed signif- 
icant light. That would be compelling evidence that Qatar-2c 
was actually a star. No such trend is apparent in the bisector 
spans. There is also no correlation with the period of Qatar- 
2b, supporting the planet interpretation. 

3.4. Stellar Parameter Classification (SPC) 

The TRES spectra were also used to classify the stellar 
parameters of Qatar-2 using SPC, a new tool (Buchhave et 
al. in preparation) for comparing an observed spectrum with 
a library of synthetic spectra. SPC has its origins in the 
procedures developed for the analysis of spectra obtained 
with the CfA Di gital Speedometers (cf. ICarnev et al.lll987l: 
Latha m et alj2002l) . It is designed to solve simultaneously for 
effective temperature (T e s), metallicity ([m/H]), surface grav- 
ity (logg), projected rotational velocity (vsini), and radial ve- 
locity (RV), taking advantage of the higher quality and more 
extensive wavelength coverage of the spectra produced by 
modern CCD echelle spectrographs. In essence, SPC cross- 
correlates an observed spectrum with a library of synthetic 
spectra for a grid of Kurucz model atmospheres and finds 
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FIG. 2. — The velocity curve for the orbit of Qatar-2b (on the bottom left, phased to the period from the orbital fit) and for the quadratic residuals due to Qatar-2c 
(on the bottom right), together with the corresponding observed velocities. The top panel shows the original velocities together with the combined velocity model. 
The panels immediately under the velocity curves show the residuals from the fits. The variations in the bisector spans are plotted in the very bottom panels. 



the stellar parameters by determining the extreme of a multi- 
dimensional surface fit to the peak correlation values from the 
grid. 

For the SPC analysis of Qatar-2 we shifted the 44 observed 
spectra to a common velocity and co-added them, to obtain a 
combined observed spectrum with SNRe of 175 in the con- 
tinuum near the Mg b features at 518 nm. An SPC anal- 
ysis yielded T efi = 4610 ±50 K, logg = 4.65 ±0.10 (cgs), 
[m/H] = -0.02 ±0.08, and vsim = 2.8 ±0.5 kms" 1 . The er- 
rors quoted here are our best guesses at the limit set by the sys- 
tematic errors, which dominate the internal errors estimated 
from the consistency of the results for the individual observa- 
tions. These values of T e ff and [m/H] are inputs to the determi- 
nation of the mass and radius of the host star (not to mention 
its age) using stellar models, as described in section 6. This 
is an important issue, because the accuracy of the mass and 
radius determined for a transiting planet is often limited by 
uncertainties in the mass and radius of its host star. 



4. FOLLOW-UP PHOTOMETRY 

4.1. KeplerCam Observations 

Light curves for four transits of Qatar-2b were obtained by 
use of KeplerCam on the 1.2 m telescope at the Fred Lawrence 
Whipple Observatory on Mount Hopkins, Arizona. Kepler- 
Cam utilizes a single Fairchild 486 4K x 4K CCD to cover 
an area of 23' x 23' on the sky, with a typical FWHM of 2"5 
for stellar images. We observed transit events of Qatar-2b on 
the nights of 2011 February 26, 2011 March 6, 2011 March 
14, and 2011 March 18. The number of images captured for 
each of these events was 79, 109, 94, and 168 respectively. 
An SDSS /-band filter was used for three of the transits, and 
an SDSS g-band filter was used for the fourth. For one of the 
/-band light curves we only covered half a transit, due to in- 
creasing cloud cover around the time of egress. The g-band 
light curve was acquired in order to look for possible color 
effects due to light from additional stars contaminating the 
Qatar-2 image, such as a background eclipsing binary. 
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TABLE 1 

Relative radial velocities of Qatar-2. 



BJD 


RV a 


TRV b 


BS C 


°"BS 




(ms- 1 ) 


(ms" 1 ) 


(ms- 1 ) 


(ms" 1 ) 



23.2 
22.6 
16.8 
17.1 
13.8 
18.6 
16.0 
21.0 
25.4 
12.3 
18.0 
25.2 
14.9 
25.2 
15.8 
23.1 
13.6 
17.1 
12.6 
20.6 
11.6 
19.1 
15.3 
19.4 
10.9 
17.5 
21.5 
15.6 
21.7 
14.6 
48.5 
36.9 
13.1 
14.9 
20.3 
19.0 
17.6 
12.5 
28.5 
15.8 
25.0 
23.6 
22.1 
28.9 



a The zero-point of these velocities is relative to the observation obtained on BJD 
2455646. The absolute velocity of that observation on the IAU system is —23.8 ± 
O.lkms -1 

k Internal errors, summed in quadrature with the uncertainty of the run-to-run zero point 

shifts, assumed to be 5 m s _1 . 

c The zero point of the bisector spans is arbitrary. 

The goal of the KeplerCam observations was to produce 
high-quality light curves, model them, and derive values for 
the radius ratio of the planet to the star, Rp/R+; the scaled 
semi-major axis of the orbit, a/R+; and orbital inclination, ; 
(or equivalently the impact parameter, b). Differential aper- 
ture photometry was carried out on the images, after the usual 
steps of bias subtraction and flat-fielding. We experimented 
with the sizes used for the aperture for the stellar images and 
the annulus used for the sky subtraction. A smaller aperture 
reduces the contribution of sky background to the noise, but 
risks larger systematic errors due to imperfect centering of the 
stellar image. We chose the aperture and annulus sizes that 
gave the best balance between these two competing sources 
of error; 8" to 9" in diameter for the stellar aperture, and an 
annulus between 30" and 60" for the sky. Sections of the out- 



TABLE 2 
Key Orbital Parameters 



Parameter Quadratic Systemic" 



P b (days) 1.337169±0.000076 1.337148 ± 0.000095 

K b (ms- 1 ) 559 ±6 559 ±6 

M b (Mi) 2.48 ±0.03 2.49 ±0.05 

e b 0.01 ±0.01 0.003 ±0.008 

/', (days ) 332:: 

AV (ins i 301:: 

\t. (Mr) 8.4:: 

<; 0.09:: 

A7 b (ms- 1 ) -282.6±5.2 

ih ,// nil- d ! 2.76 : I). 16 

d 2 j/dt 2 (ms-'d- 2 ) 0.0875±0.0046 



a This is a typical solution from the Systemic Console; we do not quote errors for the 
parameters of the outer orbit, because the Systemic Console yielded several plausible 
solutions, with differences much larger than the internal error estimates. The parameters 
given here for the outer companion should be treated with extreme caution, 
k A7 is the offset of the center-of-mass velocity for the set of relative velocities used 
for the orbital solution and reported in TablefTI 

of-transit light curves that showed contamination by incoming 
clouds or dawn were rejected. Nine reference stars were iden- 
tified, and the sum of their light curves was divided into the 
light curve for Qatar-2 in order to correct for atmospheric and 
instrumental effects. A linear trend was then fit to the out- 
of-transit sections of the resulting differential light curve and 
was used to normalize the light curve for Qatar-2b to unity. 

Figure [3] shows the KeplerCam light curves together with 
the model fits described in section 6. The residuals of the 
photometric data from the model light curves are plotted be- 
low the light curves and show very little correlated noise. In 
addition, plots of the data against column and row position on 
the CCD showed no obvious correlation, and thus are not re- 
produced here. The final photometric data for the KeplerCam 
light curves are provided in full in the electronic version of 
the journal. 

4.2. LCOGT Observations 

Four partial transit events were observed with a 0.8m tele- 
scope operated by the Las Cumbres Observatory Global Tele- 
scope (LCOGT) at the Byrne Observatory at Sedgwick Re- 
serve (BOS) near Santa Ynez, CA. The Sedgwick telescope is 
equipped with a Santa Barbara Instrument Group STL-6303E 
camera utilizing a 3K x 2K Kodak Enhanced KAF-6306E 
CCD with a pixel scale of 0.572" per pixel (2x2 binning), 
and a 14.7' x 9.8' field of view. We observed in the SDSS r 
band with exposure times of 120s. 

The images were reduced using standard routines for bias 
subtraction, dark current subtraction, and flat-field correction. 
We extracted instrumental fluxes for the stellar images using 
PyRAF and aperture photometry. Relative light curves were 
produced by dividing the fluxes of the target star by the sum 
of the fluxes of four comparison stars in each image. Each 
transit event was normalized to the out-of-transit flux on that 
particular night. Julian dates of mid-exposure were recorded 
during the observations, and later converted to BJD TDB us- 
ing the online versions of the tools described by Eastman et al. 
(2010). Aperture sizes, between 4" and 7" depending on the 
image quality, were chosen in order to minimize the photo- 
metric residuals in the resulting light curves. The final photo- 
metric data for the LCOGT light curves are provided in full in 
the electronic version of the journal. 
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FIG. 3. — KeplerCam light curves for the transits of Qatar-2b. The second light curve from the top was obtained with an SDSS g filter, the other three were 
obtained with an SDSS ;' filter. The residuals from the model are shown underneath the light curves. The exposure times were 120 seconds for the g band and 60 
for the i band. 
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5. HIGH-RESOLUTION IMAGES 

A common source of astrophysical false positives for planet 
candidates discovered by wide-angle photometric surveys is 
contamination of the target image by an eclipsing binary, ei- 
ther by a physical companion in a hierarchical system or by 
an accidental alignment with a background binary. Ground- 
based surveys that utilize cameras with short focal lengths are 
particularly vulnerable to this problem, because the detector 
pixels typically span 10" or more on the sky. The 400mm 
focal-length Canon lenses used by the QES cameras produce 
images with a typical FWHM of 7 ."5. This is roughly half the 
size of the images produced by the 200mm Canon lenses used 
by several ground-based surveys, and thus the rate of contam- 
ination by background eclipsing binaries should be four times 
better. The images of the KeplerCam follow-up photometry 
have a typical FWHM of 2."5, so this should provide nearly 
an order of magnitude additional improvement in the contam- 
ination rate. 

To push to even lower limits, we obtained high-resolution 
images with the Danish 1.54m telescope at the European 
Southern Observatory on La Silla, by use of an Andor 
iXon EM + 897 camera with an EMCCD chip, often called 
a Lucky-Imaging-camera because of its ability to obtain 
diffraction-limited images by recording lots of short expo- 
sures and collecting together the best images. The camera has 
a pixel size corresponding to 0"09, and was read at a frame 
rate of 10 Hz. In order to reach diffraction-limited images, 
one would usually stack the few percentage best quality im- 
ages from a sequence, but here we have just applied the shift 
and add technique to all the images in a sequence of 1,000 in- 
dividual exposures, which typically will reduce the seeing by 
a factor of about 3 compared to traditional CCD observations 
with the same total exposure time. 

Figure |4] shows two of the resulting images. Qatar-2 is the 
bright star in the upper right part (i.e. north-west) of the im- 
ages. The two fainter stars in the lower part of the image 
approximately 36" south-west of Qatar-2 are separated by 
slightly less than 0.6". The FWHM of the image of Qatar- 
2 is about 0."5, so these images reduce the area on the sky that 
could be contaminated by a background eclipsing binary by 
nearly another order of magnitude compared to the Kepler- 
Cam images. Although these images reduce the chance of a 
false positive due to a background eclipsing binary to a neg- 
ligible level, they do not reduce significantly the chance of a 
false positive due to contamination in a hierarchical system, 
because the angular separation of most such systems is below 
the resolution of these images. 

6. PHYSICAL CHARACTERISTICS OF THE STAR AND INNER 
PLANET 

The stellar mass and radius for Qatar-2 were initially es- 
timated using the values of T e g and [m/H] derived from the 
spectra using SPC and the value for a/R* derived from the 
KeplerCam light curves, together w ith isochrones from the 
Yonsei-Yale series of stellar mode ls (|Yi et al . 2001), follow- 
ing the procedures described by Torr es et al.l (120081) . This 
yielded a mass of M 4 = 0.742 ± 0.35M©, which was then 
imposed as an MCMC parameter controlled by a Gaussian 
prior on the simultaneous fit to the radial velocities and all 
the light curves from the QES, KeplerCam, and LCOGT pho- 
tometry. All of our exploratory fits to the radial velocities 
indicated that the orbital eccentricity for Qatar-2b was indis- 
tinguishable from circular, so in our final model we fixed that 
eccentricity to zero and fit the residuals of the radial veloci- 



ties from a circular orbit using a quadratic Taylor expansion. 
T he transit light curves w ere modeled using the formulation 
of iMandel & Agoll ((2002) in the small-planet approximation. 
A four-coefficient nonlinear limb-darkening model was used, 
employing fixed coefficients appropriate to the R band for the 
QES photometry, the SDSS r band for the LCOGT photome- 
try, and the SDSS i or g band for the KeplerCam phot ometry. 
These coefficients were determined from the tables of Claret 
(2004), interpolated to the values of T e ff and [m/H] determined 
from the TRES spectra using SPC. 

The parameter optimization was performed using 
the current version o f the MCM C code descr i bed by 
ICollier Cameron et all (120071) and iPollacco et all (120081) . 
The transit light curve is modeled in terms of the epoch 
To of mid transit, the orbital period P, the ratio of the radii 
squared d = (Rp/R+) 2 , the approximate duration tj of the 
transit from initial to final contact, and the impact parameter 
b = (acosO/R*- The radial-velocity model is defined by the 
stellar orbital velocity semi-amplitude due to the inner planet 
Kb and three coefficients in a quadratic Taylor expansion for 
the acceleration of the star due to the outer planet, the first 
coefficient being the offset of the center-of-mass velocity for 
the relative velocities used for the orbital solution, A7. The 
values of r e g- and [m/H] were treated as additional MCMC 
parameters, constrained by Gaussian priors with mean values 
and variances as determined from the TRES spectra using 
SPC. The final light curves and corresponding fits from the 
global analysis are shown in Figure [5] for the QES discovery 
photometry and for the KeplerCam and LCOGT follow-up 
photometry. The correlations between the various posterior 
parameter estimates are shown in the matrix of plots in 
Figure [6] The final physical and orbital parameters and error 
estimates for the star and planet are reported in Table [3] It 
is reassuring to see the good agreement between the stellar 
parameters estimated from the spectra using SPC and the 
final values from the global analysis: T e s = 4650 ± 50 K from 
SPC compared to 4645 ± 50, and logg = 4.65 ± O.lOcgs 
versus 4.601 ±0.018. 

7. DISCUSSION 
7.1. Hot Jupiters with Companions 

One of the most important results from the Kepler mission 
is the discovery of a rich population of multiple transiting 
planets in sy stems whose orbits must be coplanar with in a de- 
gree or two (lLatham et al.ll201 ll iLissauer et al.ll201 hV Espe- 
cially striking is the very low occurrence rate for hot Jupiters 
in the multiple systems found by Kepl er. In the sample o f 
1235 planet candidates announced by iBorucki et al.l d201 ll) . 
there are 117 systems that contain a planet candidate more 
massive than Saturn (Mp > 0.3 Mj) in an orbit with period 
shorter than 10 days, but only five of these systems harbor 
an additional transiting planet candidate. The remaining 112 
transiting hot Jupiters are all in singles. 

The occurrence rate is also very low for companions to 
hot Jupiters found by ground-based surveys, both photometric 
and spectr oscopic. A revi ew of the Extrasolar Planet Ency- 
clopaedia ( Schneiderl l201 ll) reveals that Qatar-2b joins HAT- 
P-13b, HAT-P-17b, and HAT-P-31b as the fourth transiting 
hot Jupiter with a confirmed outer planet, while only four 
hot Jupiters found by radial velocities have outer compan- 
ions (HD 217107b, HIP 14810b, HD 187123b, and v Andb). 
Of course, additional radial-velocity monitoring may reveal 
other hot Jupiters with outer companions, and there are al- 
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FIG. 4. — Lucky camera images of Qatar-2. North is up and west is left. The faint pair of stars 36" to the south-west are separated by slightly less than 0"6 



TABLE 3 

Planetary and Stellar Parameters. 



Parameter 


Symbol 


Value 


Unit 


Transit epoch (BJD TDB) 


To 


2455624.26679 ±0.00011 


days 


Orbital period 


1\ 


1.3371 182 ±0.0000037 


days 


Transit duration 


t T 


0.07540 ±0.00049 


days 


Planet/star area ratio 


(Rp/R*) 2 


0.02725 ± 0.00040 




Impact parameter 


b 


0.19±0.10 


R* 


Scaled stellar radius 


R*/a 


0.1541 ±0.0030 




Stellar density 


P* 


2.05 ±0.12 


PO 


Stellar effective temperature 


T c ff 


4645 ±50 


K 


Spectroscopic Metallicity 


[m/H] 


(fixed) 




Stellar surface gravity 


logg* 


4.601 ±0.018 


(cgs) 


Projected stellar rotation speed 


v sin i 


2.8 (fixed) 


kms - 


Stellar radius 


R* 


0.713±0.018 


Re 


Stellar mass 


M, 


0.740 ±0.037 


M e 


Orbital separation 


a 


0.02149 ±0.00036 


AU 


Orbital inclination 


i 


88.30 ±0.94 





Stellar reflex velocity 




558.7 ±5.9 


ms 


Center-of-mass velocity offset 


A7 


-282.6 ±4.0 


ms~' 


Drift in center-of-mass velocity 


dj/dt 


-2.74 ±0.15 


ms 1 


Quadratic velocity drift term 


d 2 -y/dt 2 


0.0870 ±0.0043 


ms~'c 


Orbital eccentricity 


e 


(fixed) 




Planet radius 


Kp 


1.144 ±0.035 


Ri 


Planet mass 


M p 


2.487 ±0.086 


M s 


Planet surface gravity 


l°g£p 


3.638 ±0.022 


(cgs) 


Planet density 


Pp 


1.66 ±0.13 


PJ 


Planetary equilibrium temperature 


Tv 


1292 ±19 


K 



ready hints of velocity drifts for WASP-8 and WASP-22. All 
of the confirmed outer companions have relatively long or- 
bital periods (see Table|4j>, so it is not surprising that none of 
them have shown transits, so far. Nevertheless, it is tempting 
to speculate that a close-in giant planet may stir up the orbits 
of other inner planets in its system, while a system of small 
planets is more likely to preserve the flatness of the disk from 
which it formed, allowing small planets to survive in surpris- 
ingly compact configurations. 

In only one of the systems (HIP 14810) listed in Table|4]is 
the outer companion less massive than the hot Jupiter. In the 
rest of the systems the nearest companions are more massive, 
often by large factors. 



7.2. Orbits ofQatar-2b andHAT-P-31b 

Why i s the orbit of HAT-P- 3 lb eccentric, with e = 0.2450 ± 
0.0045 (Kippin g et al J 120 111) , while the orbit of Qatar-2b is 
indistinguishable from circular? One possible explanation is 
that both orbits started out with significant eccentricity, per- 
haps as the result of a dynamical encounter that sent each 
planet on a path close to its host star, and the orbit of Qatar-2b 
has since been circularized by tidal forces, while the orbit of 
HAT-P-31bhasnot. 

The circularization time scale goes something like 
(Mp/M*)(a/R±) 5 . HAT-P-31b has a period of P = 5.00 days 
and mass of Mp = 2.17 'Mj, orbits a late F star with mass 
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FIG. 5. — The fit of the light curves found by the global analysis. The first curve in the upper panel shows the QES discovery light curve in the R band; 
the second shows the KeplerCam SDSS g-band light curve; the third shows the LCOGT SDSS r-band light curve, and the fourth shows the KeplerCam SDSS 
i-band light curve. The lower panel shows the corresponding residuals from the global fit. The dip in the QES light curve may be too shallow due to aggressive 
detrending. 



M* = 1.22M©, a/R+ = 8.9, and is about 3 Gyr old. In con- 
trast, Qatar-2b has a much shorter period of Pb = 1-34 days, 
similar mass of Mp = 2.49 Mj, orbits a K dwarf with mass of 
= 0.74M Q , a/R+ = 6.3, and is very likely a member of 
the Galactic disk. If Qatar-2 is three times older than HAT-P- 
31, the system age divided by the circularization time scale 
is nearly an order of magnitude longer for Qatar-2b, long 
enough so that there has been time for tidal forces to circu- 
larize the orbit. In addition, the host star for HAT-P-31b is 
close to the mass where the outer envelope no longer supports 
a convection zone, and the circularization time scale is even 
longer than implied by the simple relation that we adopted for 
the sake of this discussion. 



Another effect that needs to be explored more carefully 
with dynamical simulations is the amount of eccentricity that 
Qatar-2c can pump into the orbit of Qatar-2b. Are the tidal 
forces working to circularize the orbit of Qatar-2b strong 
enough to keep the orbit circularized despite the perturbations 
from Qatar-2c? Finally, it would be interesting to estimate the 
size and patterns of the transit time variations expected for 
Qatar-2b due to perturbations by Qatar-2c. 

7.3. The Value of Deep Transits 

With a depth of about 3.5%, the transits of Qatar-2b are 
deeper than the transits of any other planet listed in The Ex- 
trasolar Planets Encyclopaedia (Schneider 2011). The clos- 
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FIG. 6. — The correlations between the various posterior estimates from the MCMC global fit. 
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TABLE 4 
Hot Jupiters with Companions 



Planet 


M v (M,r 




P(d) 


Qatar-2b 


2.49 ±0.09 


1.14 ±0.04 


1.337118 


c 


8.4:: 




332:: 


HAT-P-13b 


0.85 ±0.04 


1.28 ±0.08 


2.916243 


c 


14.5 ±1.0 




448.2 ±1.0 


HAT-P-17b 


0.53 ±0.02 


1.01 ±0.03 


10.338523 


c 


1 4 +L1 




1798+| 


HAT-P-31b 


2.2±0.1 


1.1 ±0.4 


5.005424 


c 


>3.4 




> 1022 


HD 217107b 


1.85±0.05 




7.12689 


c 


2.5 ±0.3 




4210 ±190 


HIP 14810b 


3.9 ±0.3 




6.673855 


c 


1.3 ±0.1 




147.73 


d 


0.57 ±0.05 




962 ±15 


HD 187123b 


0.52 ±0.04 




3.0965828 


c 


2.0 ±0.3 




3810±420 


v Andb 


0.69 ±0.03 




4.617136 


c 


10.3+ff 
1.06 ±0.03 




240.94 


d 




1282 


e 




3848 



a (AU) 


e 




M,(M Q ) 


0.0215 
0.82:: 
0.0426 

1.19 
0.088 

2.8 
0.055 


(fixed) 
0.09:: 
0.014 ±0.05 
0.67 ±0.02 
0.346 ±0.007 

1+ 02 
0.245 ±0.005 


88.30 ±0.94 
83.3 ±0.6 
89.2 ±0.2 

87.12? 


0.74 
1.22 
0.86 

1.22 


0.073 
5.27 
0.069 
0.55 

1.9 
0.043 

4.9 
0.059 
0.83 
2.53 
5.25 


0.132 ±0.005 
0.52 ±0.03 
0.1427 ±0.0009 
0.16±0.01 
0.17 ±0.04 
0.01 ±0.01 
0.25 ±0.03 
0.013±0.016 
0.245 ±0.006 
0.316±0.006 
0.0054 ±0.0004 


>30 
8±1 
24±1 


1.02 
0.99 

1.06 
1.27 



a For the transiting hot Jupiters (the first four entries) the actual mass is reported; for the hot Jupiters discovered by Doppler surveys (the final four entries), the minimum mass is listed, 
except for v And c&d, which are actual masses enabled by inclinations determined with the Fine Guidance Sensors on HST. 
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est rival is CoRoT- 2b, with reported transit depths of 3.2% 
dAlonso et alj|2008l) . CoRoT-2 is a late G dwarf with T eB = 
5625 K, while Qatar-2 is a late K dwarf with r eff = 4645 
K; the smaller radius of Qatar-2 is the main reason for the 
deeper transits. A more interesting comparison may be the 
bright (V = 7.7 mag) early K dwarf HP 189733 wi th T eB = 
5050 K (cf. iBouchv et all 12001 iBdcoseFall l2006h . which 
has been a f avorite target for studies of the planet's at- 
mosphere ("cf . iKnutson et all 120071: ISwain. Vasisht. & Tinettil 
120081: iGrillmair et al.1 120081: iGibson. Pont & Aigrainl 120 iW . 
The contrast between the planet and host star is even more 
favorable for Qatar-2, but the system is much fainter (V = 
13.3 mag). Thus HD 189733 is likely to continue as a 
top target for follow-up studies of hot Jupiters, although 
it may eventually be joined by other nearby bright transit- 
ing systems discovered by all-sky transit surveys such as 
TESS or PLATO, or by targeted searches of small cool 



stars such as MEarth. Nevertheless, Qatar-2 is a good tar- 
get for amateurs and outreach projects such as MicroObser- 
vatory (http://mo-www.harvard.edu/MicroObservatory/) with 
telescopes of modest size but covering wide fields of view, 
because of the availability of many more reference stars of 
comparable magnitude. 
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